Feedback-stabilization of an arbitrary pure state of a two-level atom 
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Unit-efficiency homodyne detection of tlie resonance fluorescence of a two-level atom collapses 
the quantum state of the atom to a stochastically moving point on the Bloch sphere. Recently, 
Hofmann, Mahler, and Hess [Phys. Rev. A 57, 4877 (1998)] showed that by making part of the 
coherent driving proportional to the homodyne photocurrent can stabilize the state to any point on 
the bottom half of the sphere. Here we reanalyze their proposal using the technique of stochastic 
master equations, allowing their results to be generalized in two ways. First, we show that any 
point on the upper or lower half, but not the equator, of the sphere may be stabilized. Second, we 
consider non-unit-efficiency detection, and quantify the effectiveness of the feedback by calculating 
the maximal purity obtainable in any particular direction in Bloch space. 
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I. INTRODUCTION 

Although classical models of feedback schemes have 
been used for a long time to control dynamical noise, an 
analogous quantum theory of feedback has been devel- 
oped only in the last fifteen years Recently there 
has been considerable interest in quantum feedback as a 
way to fight decoherence in isolated quantum systems, 
using the approach of Refs. The central idea is to 
use a continuous measurement record, whose existence is 
due to the couphng of the system to a bath, to control 
the dynamics of the system so as to counteract the noise 
introduced by that bath and possibly other baths. For 
example, it has been suggested as a way to create optical 
squeezed states 1^, to create micromaser number states 

, to correct errors in quantum bits jl^] , and to protect 
optical and microwave Schrodinger cat states against dis- 
sipation (ll]-[l3). 

Decoherence in quantum systems can be loosely de- 
fined as loss of purity. Therefore the ultimate success in 
using feedback to fight decoherence would be to create an 
arbitrary stable pure state in the presence of dissipation. 
This goal was realized (better even that they realized) 
by Hoffman, Mahler and Hess (HMH) for a very 

simple system: a resonantly driven two-level atom. They 
showed that by using the photocurrent derived from unit- 
efficiency homodyne detection of the atom's fluorescence 
to control part of the driving field of the atom, it is possi- 
ble to exactly cancel the noise introduced by the electro- 
magnetic vacuum field when the atom is in a particular 
pure state. By choosing the driving strength and feed- 
back strength appropriately, any pure state on the Bloch 
sphere may be picked out, although HMH claimed that 
only pure states on the lower half of the sphere would be 
stable under their scheme. 

HMH chose to describe detection and feedback in their 
system in a way different from (but equivalent to) the 
standard approach in Refs. In this paper we refor- 

mulate their theory using the latter approach. This has 



the advantage of enabling a number of generalizations of 
their results. First, we revisit the question of stability 
and find that, contrary to the claims of HMH, the states 
in the upper half of the Bloch sphere can be stabilized 
as well as those in the lower half (this is what was bet- 
ter than they realized). The only states which cannot 
be stabilized, in the sense that an arbitrary initial state 
would not always end up in the desired state, are those 
on the equator of the Bloch sphere; that is, those which 
are equal superpositions of excited and ground states. 

Our second generalization is to consider how effective 
feedback is with 77 < 1; that is, with non-unit-efficiency 
detection. In this case it is not possible to stabilize the 
atom at any pure state, except the ground state which 
is trivially stable by setting the driving and feedback to 
zero. Instead, we aim to produce a steady state which 
is as close as possible to a given pure state. For the 
two-level atom, this is equivalent to trying to create a 
state which is as pure as possible in a particular direc- 
tion in Bloch space. Not surprisingly (given the above 
result), we find that states near the equator cannot be 
well-protected against decoherence. We also find an echo 
of the distinction HMH found between the upper and 
lower halves of the Bloch sphere, in that states in the 
upper half sphere are affected much more by loss of de- 
tection efficiency that those in the lower half. 

The paper is organized as follows. In Sec. II we 
present the model of a driven two-level atom, includ- 
ing the stochastic Schrodinger equation for unit-efficiency 
homodyne detection. In Sec. Ill we use this equation to 
derive the driving and feedback required to stabilize the 
atom in an arbitrary pure state. These results agree with 
those of HMH. However, our stability analysis disagrees 
substantially with theirs. In Sec. IV we present entirely 
new analytical results relating to the effect of non-unit- 
efficiency detection. In Sec. V we give numerical simu- 
lations of the stochastic evolution equations, illustrating 
the issues discussed in the preceding two sections. In 
Sec. VI we summarize and interpret our results, explain 
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their significance, and discuss the possibihty of future 
work. 

II. HOMODYNE DETECTION 
A. Master Equation 



FIG. 1. Diagram of the experimental apparatus. The laser 
beam is split to produce both the local oscillator (3 and the 
field ao which is modulated using the current I{t). The mod- 
ulated beam, with amplitude proportional to a-|-A7(t), drives 
an atom at the centre of the parabolic mirror. The fluores- 
cence thus collected is subject to homodyne detection using 
the local oscillator, and gives rise to the homodyne photocur- 
rent I{t). 



Consider a atom, with two relevant levels {jg), |e)} and 
lowering operator a = |5)(e|. Let the the decay rate be 
7, and let it be driven by a resonant classical driving 
field with Rabi frequency 2a. This is as shown in Fig. ||, 
where for the moment we are omitting feedback by set- 
ting A = 0. This system is well- approximated by the 
master equation 



(2.1) 



It is easy to sh ow that the stationary solution of the 
is 
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where the Lindblad |lq | superoperator is defined as usual 
V[A]B = ABA^ - {A^A,B}/2. In this master equation 
we have chosen to define the — a+a^ and ay = ia—ia^ 
quadratures of the atomic dipole relative to the driving 
field. The effect of driving is to rotate the atom in Bloch 
space around the y-axis. The state of the atom in Bloch 
space is described by the three-vector {x,y,z). It is re- 
lated to the state matrix p by 
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(/ + xax + ycTy + zffz) . 



(2.2) 



For 7 fixed, this is a family of solutions parameterized by 
the driving strength a G (— cx),cx)). All members of the 
family are in the x-z plane on the Bloch sphere. Thus 
for this purpose we can reparametrize the relevant states 
using r and 9 by 



X = r sm fc 
z — r cos f 



(2.6) 
(2.7) 




where G [— tt, tt]. Since 



Tr[p'] = - (1 + x2 + ,/ + z2) (2., 



is a measure of the purity of the Bloch sphere, r — 
\/a;2 + z2j the distance from the centre of the sphere, is 
also a measure of purity. Pure states correspond to r = 1 
and maximally mixed states to r = 0. 

The locus of solutions in this plane (an ellipse) is shown 
in Fig. ^. Since Zgs < 0, all solutions are in the lower 
half of the Bloch sphere. That is, we are restricted to 
16*1 > 7t/2. Also, it is evident that the smaller \0\ is (that 
is, the more excited the atom is), the smaller r is (that 
is, the less pure the atom is). At \9\ = tt the stationary 
state is pure, but this is not surprising as it is simply the 
ground state of the atom with no driving. As \6\ Tr/2 
we have r — > 0. This can only be approached asymptoti- 
cally as |a| — > c». In summary, the stationary states we 
can reach by driving the atom are limited, and generally 
far from pure. 
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FIG. 2. Locus of the solutions to the Bloch equations. The 
ellipse in the lower half plane is the locus for the equations 
with driving only. The full circle (minus the points on the 
equator) is the locus for the equations with optimal driving 
and feedback, as defined in Sec. III. 



B. Homodyne Detection 

Now consider subjecting the atom to homodyne de- 
tection. As shown in Fig. we assume that all of the 
fluorescence of the atom is collected and turned into a 
beam (represented in Fig. |^ by placing the atom at the 
focus of a mirror). Ignoring the vacuum fluctuations in 
the field, the annihilation operator for this beam is ^/7cr, 
normalized so that the mean intensity 7 ^cr^cr) is equal 
to the number of photons per unit time in the beam. 
This beam then enters one port of a 50:50 beam split- 
ter, while a strong local oscillator /? enters the other. To 
ensure that this local oscillator has a fixed phase relation- 
ship with the driving laser used in the measurement, it 
would be natural to utilize the same coherent light field 
source in the driving process and as the local oscillator 
in the homodyne detection. This is as shown in Fig. |l|. 

Again ignoring vacuum fluctuations, the two field op- 
erators exiting the beam splitter, h\ and 62, are 



fefe = Ul^ - (-1)'/?] /%/2. 



(2.9) 



When these two fields are detected, the two photocur- 
rents produced have means 

4 = (|/3P - (-1)' (V7/3'^^ + V7^/3*) + 10'^ o) /2. 

(2.10) 

The middle two terms represent the interference between 
the system and the local oscillator. 

The ideal limit of homodyne detection is when the lo- 
cal oscillator amplitude goes to infinity, which in prac- 
tical terms means |/3p ^ 7. In this limit, the rate of 



the photodetections goes to infinity and thus it should 
be possible to change the point process of photocounts 
into a continuous photocurrent with white noise. Also, 
the only relevant quantity is the difference between the 
two photocurrents. Suitably normalized, this is ]l7| , |l8t 



hit) -hit) 



(2.11) 



A number of aspects of Eq. ( ^.11 ) need to be explained. 
First, $ — arg/3, the phase of the local oscillator (de- 
fined relative to the driving field). Second, the subscript 
c means conditioned and refers to the fact that if one 
is making a homodyne measurement then this yields in- 
formation about the system. Hence, any system aver- 
ages will be conditioned on the previous photocurrent 
record. Third, the final term ^(t) represents Gaussian 
white noise, so that 

^it)dt = dWit), (2.12) 

an infinitesimal Wiener increment defined by | |T9|] 



dWit)^ = dt, 
E[dW^(t)] = 0. 



(2.13) 
(2.14) 



Since the stationary solution of the master equation 
confines the state to the x-z plane, it makes sense to 
follow HMH by setting <I> = 0. In that case. 



(2.15) 



That is, the deterministic part of the homodyne pho- 
tocurrent is proportional to Xc = {<^x)c- This should be 
useful for controlling the dynamics of the state in the 
x-z plane by feedback, as we will consider in Sec. III. 
Of course, all that really matters here is the relationship 
between the driving phase and the local oscillator phase, 
not the absolute phase of either. 

The conditioning process referred to above can be 
made explicit by calculating how the system state 
changes in response to the measured photocurrent. As- 
suming that the state at some point in time is pure (which 
will tend to happen because of the conditioning anyway), 
its future evolution can be described by the stochastic 
Schrodinger equation (SSE) 0,|l8| 

d\Mt)) = A,it)\Mt))dt + B,it)\^,it))dWit). (2.16) 

This is an Ito stochastic equation |^ with a drift term 
and a diffusion term. The operator for the drift term is 



Mt) = I [-aV + (a,)^ {t)a - (t)/4 



while that for the diffusion is 



(2.17) 
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(2.18) 



Both of these operators are conditioned in that they de- 
pend on the system average 



{cT,),{t) = {Mt)\^x\Mt))- 



(2.19) 



O n average, the system still obeys the master equation 
( ^.l|) . This is easiest to see from the stochastic master 
equation (SME), which allows for impure initial condi- 
tions. The SME can be derived from the SSE by con- 
structing 

d{\A)m) =(d|^c))(V'c| + |V^c)(rf(V'c|) 

-t- (d|V^c))(d(V^c|), (2.20) 



using the Ito rule ( 2.13| ), and then identifying \ipc){'^pc\ 
with pc- The result is 

dpc = dt'^'D[ij]pc - idta[ay,pc] + dW{t)^TL[(7]pc, 

(2.21) 

where n[A\B = AB + BA^ - Tt[AB + BA^. Although 
this has been derived assuming pure initial conditions, 
it is valid for any initital conditions [Q. This is also an 
Ito equation, which means the evolution for the ensemble 
average state matrix 



p{t) = E[pc(t)] 



(2.22) 



is found simply by avera ging over the photocurrent noise 
term by using Eq. ( 2.14). This procedure yields the orig- 
inal master equation ( ^.l[ ) again. The general term for 
the stochastic conditioned evolution of the system, be it 
described by a SSE or SME, is a quantum trajectory p7| , 
and the quantum trajectory is said to unravel the master 
equation |17|. 



III. FEEDBACK WITH UNIT-EFFICIENCY 
DETECTION 

A. SSE including feedback 



the form (2.16). The effect of the feedback Hamiltonian 
can be shown to change the drift and diffusion op- 
erators to 

+ AV7 {'^.)c {t)^y - 2'^^^] - ^V2, (3.2) 



Beit) =^[(J- {(Jx), {t)/2] - iXay 



(3.3) 



Say we wish to stabilize the pure state with Bloch an- 
gle 9, as defined in Eqs. ( |2.6| ) and (2/7), with r = 1 of 
course. In terms of the ground and excited states, this 
state is 



cos-|e) -|-sin-|5). 



(3.4) 



Now for this state to be stabilized we must have 

Ac{t)dt + ]3c{t)dW{t) \ \e) oc 16*). (3.5) 

We cannot say the left-hand-side should equal zero be- 
cause a change in the overall phase still leaves the phys- 
ical state unchanged. However, we can work with this 
equation, and simplify it by dropping all terms propor- 
tional to the identity operator in Ac{t) and Bc{t). We 
can also demand that it be satisfied for the deterministic 
and noise terms separately, because dW(t) can take any 
value. This gives the two equations 



[7 (^^a^cr + sinOa^ — i2aay 
+ ^-^7 i sin 6*0-2^/2 — cr^)] 



(3.6) 
(3.7) 



where we have put {cx)^ (t) equal to sin0, its value for 
the state \6). 

Solving the first equation easily yields the condition 



A = -^(1 -(-cosff). 



(3.8) 



We now include feedback onto the amplitude of the 
driving on the atom, proportional to the homodyne pho- 
tocurrent, as done by HMH. This is as shown in Fig. |^, 
where the driving field passes through an electro-optic 
amplitude modulator controlled by the photocurrent, 
yielding a field proportional to a -I- XI (t). This means 
that the feedback can be described by the Hamiltonian 



77fb - Xayl{t). 



(3.1) 



In this paper we are assuming instantaneous feedback, 
where the time delay in the feedback loop is negligible. 

Since the homodyne photocurrent (2.11) is defined in 
terms of system averages and the noise dW{t), the SSE 
including feedback can still be written as an equation of 



This is equivalent to the feedback condition derived by 
HMH, stated as Eq. (35) of Ref. ||l|l. Substituting this 
into the second equation gives, after some trigonometric 
manipulation, the second condition 



7 

a = — sin 9 cos ( 
4 



(3.9) 



Again, this agrees with the driving strength of HMH, 
given as Eq. (44) of Ref. It is worth emphasizing 

that the derivation given here is entirely different in de- 
tail from that of HMH, and so is an independent verifica- 
tion of their result. These functions are plotted in Fig. ||. 
Note that there are two points with the same values of 
both A and a, a.t 9 = ±7r/2. 
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to be norm- preserving). It is not difficult to find these 
eigenvalues, and in terms of they are 
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FIG. 3. Plot of the optimal driving (a, solid) and feedback 
(A, dashed) required to produce a pure state with Bloch angle 
9. For this plot we have set 7 = 1 so that ct and A are dimen- 
sionless. The purity (r^, starred) is one except for 9 = ±7r/2, 
where the feedback is not stable. 



B. Stability 

The preceding derivation seems to show that any pure 
state can be stabilized by a suitable choice of driving 
and feedback. Indeed our derivation proves that that if 
one prepares a state in exactly the pure state one de- 
sires, then the feedback scheme of HMH which we have 
analyzed will keep the system in that state. However, 
to discuss stability we need to know what will happen 
for states which are not initially in the desired state. To 
deal with this it is much more convenient to use the SME 
rather than the SSE, as will become apparent. 

The SME can be constructed from the SSE in the same 
way as before. The result is 

dpc = dt-f'D[a]pc - idta[ay,pc] 

- idtX[ay, apc + Pccr^] + dt{X'^ /^)'D[<Ty]pc 



-dW{t)H[^(j-iXay]p,. 



(3.10) 



Also as before, this is an Ito stochastic equation, which 
means that the ensemble average can be found simply by 
dropping the stochastic terms. This time, the result is 
not the original master equation, but rather the feedback- 
modified master equation 



P 



-i[aay,p] +'D[^C7 ~ iXcry]p = Cp. 



(3.11) 



Here we have put the Liouvillian superoperator £ in a 
manifestly Lindblad form. 

Now we have shown already that the pure state p = 
\d){6\ must be a solution of th is maste r eq uation, for the 
appropriate choices of A (3.8) and a (3.9). But for it to 
be a stable solution we require all of the eigenvalues of 
the resulting Cg to have a negative real part (except for 
the one eigenvalue that is always zero, as required for Cg 



7/2, -7/2,-7COS^ 



(3.12) 



Evidently the state \d) will be stable for all 6 except 
6 = ±7r/2. That is, all states are stable except those 
on the eq uator. This is contrary to the conclusion of 
HMH [g5|, based on a linearized stability analysis, that 
"long term stability of . . . inverted states [i.e. states in 
the upper half plane] cannot be achieved." We emphasize 
that our stability analysis contains no approximations. 

In hindsight, the lack of stability for pure states on 
the eq uato r coul d h ave been predicted from the expres- 
sions ( |3.8|) and (3^). As discussed above and shown in 
Fig. |3[ the values for driving and feedback for 9 — n/2 are 
the same as those for 6 = —tt/2. This means that both 
p = |7r/2)(7r/2| and p = |— 7r/2)(— 7r/2| are solutions of 
Cgp — ioi 9 — tt/2 or — 7r/2. By linearity, any mixture 
of |7r/2)(7r/2| and |— 7r/2) (— 7r/2| will be a solution also. 
Hence any deviation away from one of these pure state 
will not necessarily be suppressed, and the system lacks 
stability. With random external perturbations, the sys- 
tem will eventually reach an equal mixture of |7r/2)(7r/2| 
and |— 7r/2)(— 7r/2|, which is a state with r = (minimum 
purity). This is why we have plotted a value of r = in 
Fig. ^ for \0\ = Tr/2. We also plot r as a function of 
in Bloch space in Fig H, giving the locus of states which 
can be stabilized by feedback. This can compared to the 
locus of the mixed states which are accessible by driving 
alone. We will return to the stability issue in the context 
of stochastic dynamics in Sec. V. 



IV. FEEDBACK WITH NON-UNIT-EFFICIENCY 
DETECTION 

We have seen that the stochastic master equation is a 
very useful representation of a quantum trajectory, as it 
allows the unconditioned (deterministic) master equation 
to be derived easily, and this latter equation is all that 
is required for a completely rigorous stability analysis. 
The SME is also superior to the SSE in that allows in- 
efficient detection to be described. In a real experiment 
this has to be taken into account. The effect of non-unit 
rj on feedback in the present system is of interest both 
in itself, and because of the extreme nonlinearity of the 
system dynamics as compared to other quantum optical 
feedback systems such as considered in Ref. |8|. 

As explained in Ref. |l8j, the homodyne photocurrent 
from a detection scheme with efficiency rj is 



(4.1) 



Here we have used a normalization such that the deter- 
ministic part does not depend on rj. The effect of de- 
creased efficiency is increased noise. This means that 
we can retain the same feedback Hamiltonian as above 
31), without changing the significance of the feedback 
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parameter A. The SME with 7] < 1, including feedback, 



dpc = dt-/'D[a]pc ~ idta[(7y, pc] 

- idtX[ay,apc + Pccr^] + dt{X^ /^vi)'D[ay\pc 



■ dW{t)n[y^a - iX(jy/y/Tj]p^ 



(4.2) 



The no- feedback SME, analogous to Eq. (2.21 



), can be 

obtained simply by setting A 0, and was derived in 
Ref. [HI. 

Once again, it is easiest for the moment to just con- 
sider the ensemble average evolution by averaging dW to 
zero. The Lindblad form of the resulting master equation 
is 

p = -i[aay,p] + Vly^a - iX(Ty]p + {X^ / r])V[cFy\p. (4.3) 

We do not know a priori what values of A and a to 
choose to give the best resu lts wi th inefficient detection, 
as the SSE analysis in Sec. Ill A obviously does not ap- 
ply. Hence we simply solve for the stationary matrix in 
terms of a and A. Using the Bloch representation we find 



^ss — 
Vss = 
^ss — 

where 



4ar7^(7 
0, 



2^X)/D, 

2A)(77/-f 4/777A 



(4.4) 
(4.5) 

AX^)/D, (4.6) 



D = -/"^rf -f Q-i^'^rfX + 2777(8 
A^). 



4?7)A2 



16V7^A3 

(4.7) 



The question now arises, what do we mean by "best 
results" for the feedback system. We cannot hope any- 
more to produce stable pure states anywhere on the Bloch 
sphere. However, we can pick a direction 9 on the Bloch 
sphere and ask how close can we ge t to a pure st ate? 
That is, we use the radius r in Eq. (2.6) and Eq. (2.7) 
as the quantity to be maximized, for each 9. From these 
two equations we have 



tan( 



(4.8) 



^From Eqs. (4.4) and (4.6) we can immediately find the 
desired driving in terms of A and 9 as 



-V^viV^ + 2A)(777 + 4^7/A + 4A2) tan6l 
4r;2(7 + 2V7A) ' 



(4.9) 



The aim is then, for each 9, to find the feedback A which 



maximizes 



(4.10) 



This was done numerically using matlab. 

The results of our calculations are shown in Fig. ^, 
where we plot the locus in Bloch space of the best (most 
pure) stationary states which can be achieved by feed- 
back from non-unit-efficiency detection. We use a vari- 
ety of values of rj. A number of points are worth noting. 



First, and most obviously, the degree of purity (measured 
by the r, the distance from the origin) decreases with 77. 
Second, the gap at the equator for rj — I quickly widens, 
so that the purity of the best states with 9 close to 7r/2 
is small. Third, the purity of the best states in the upper 
half of the Bloch sphere is affected much more by loss of 
detection efficiency than those in the lower half. Fourth, 
in the limit 77 = 0, the best solutions correspond to the 
no-feedback solutions shown in Fig. ^ This is not sur- 
prising, since with 77 = the photocurrent contains no 
information about the system (as the noise is infinitely 
large) and hence there is no point doing feedback. Since 
the stationary states with no feedback are confined to the 
lower half of the Bloch sphere, this explains why the best 
states with feedback in the lower half are less affected as 
r] decreases than those in the upper half. 
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FIG. 4. Locus of the solutions to the Bloch equations with 
optimal feedback for different values of detector efficiency rj. 
from the outside in, we have rj = 1, 0.8, 0.6, 0.4, 0.2, 0. 

For the particular value 77 — 0.8 we plot in Fig. || the 
values of a and A (as well as purity, quantified as r^) 
versus 9. By comparing this plot with Fig. |3| one obtains 
some idea of the effect of inefficient detection. A number 
of features remain the same. First, a is antisymmetric 
in 9, while A is symmetric. Recalling that the determin- 
istic part of the feedback is proportional to A (cx)^, the 
feedback itself is actually antisymmetric as well as the 
driving. Second, the magnitude of the feedback is zero 
for \9\ — Tx (the ground state) and increases monoton- 
ically to a maximum of ^ 9 — Q (in the direction 
of the excited state). Third, the driving is zero at the 
ground state and at = 0, and also changes sign as one 
passes through the equatorial place. The most obvious 
difference between the parameters for 77 = 1 and those for 
r\ = 0.8 parameter is that the latter have a discontinuity 
at \9\ — 7^/2. The feedback parameter A jumps as one 
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crosses the equatorial plane, while the driving a asymp- 
totes to +00 on one side and —00 on the other. These 
extreme variations in the driving do not prevent the best 
purity from approaching zero in the equatorial plane. 
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FIG. 5. Plot of the optimal driving (a, solid) and feedback 
(A, dashed) required to produce the most pure state with 
Bloch angle Q. For this plot we have set 7 = 1 so that a and 



A are dimensionless. 
starred) . 



The purity obtain is also plotted (r 



V. STOCHASTIC DYNAMICS 



A. Stochastic Bloch Equations 



(5.2) 

where k — A^/?7 + If we ignore the final (noise) 

term, we get the Bloch equations from the master equa- 
tion dJ). 



B. unit-Efficiency 

In the case r\ considered in Sec. Ill, both the de- 
terministic and stochastic dynamics disappeared in the 
steady state for the appropriate choice of a and A. Be- 
cause the stationary solution of the SSE was a unique 
pure state, that was necessarily also the stationary solu- 
tion of the master equation found by averaging over the 
noise in the equivalent SME. Thus there was no distinc- 
tion between the unconditioned and conditioned states. 
There are two exceptions to this lack of distinction. The 
first is in the transients, before the system reaches its 
steady state. The second is for the special case \B\ = 7r/2. 
In this subsection we investigate these exceptions. 
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So far we have considered the stochastic conditioned 
dynamics for the system state in order to derive the pa- 
rameters A and a such that for 77 = 1 those dynamics 
are banished in the steady state. In this section we will 
consider them in more detail, highlighting the difference 
between the r\ — \ case and the 77 < 1 case, and also 
looking in more detail at the special case of \d\ = 7r/2. 
The most convenient way to treat the stochastic dynam- 
ics in general is through the stochastic Bloch equations 
(SBE). These are simply the stochastic equations for the 
conditioned Bloch vector, defined by 
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FIG. 6. Typical quantum trajectories for optimal feedback 
with rj — 1, shown by Xc (plusses) and Zc (dots) as functions 
of time, (a) shows a single trajectory for 9 = n/6, and (b) 
two trajectories for = n/2. 

First, the transient behaviour was simulated using the 
SBE with 77 = 1. We chose the initial state to be the 
ground state, and evolved the system stochastically from 
t = to t = 107"^. With this choice of initial condition, 
Uc — for all time. We verified that in each trajectory 

+ = 1 to a good approximation (indicating a pure 
state), but that the ensemble averages over many trajec- 
tories 
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E[Ze] 



(5.3) 



obey the deterministic Bloch equations. A typical trajec- 
tory for 9 = 7r/6 is shown in Fig. ||(a). We see that the 
initial evolution is very erratic, but that on a time scale 
of a few the system relaxes towards a steady state 
which is pure and stationary. By t = 107~^ the system 
is locked in a stable pure state for all intents. We have 
also illustrated another typical trajectory in Bloch space, 
in Fig. 0. 



N 




FIG. 7. Typical quantum trajectories in Bloch space for 
t G [0, 107~^] under optimal feedback for 6 — 7r/6, starting at 
the ground state. The plusses are for 77 = 1 and the dots for 
77 = 0.8. The locus for the deterministic stationary states for 
77 = 0.8 are also shown; the relevant state for this quantum 
trajectory is at the intersection of the locus and the ray at 
6 = 7r/6. Note that the quantum trajectory for 77 = 0.8 wan- 
ders around this average position, while that for 77 = 1 stops 
precisely at the desired pure steady state. 

It is easy to verify that by putting 77 = 1 and 

(xcj/c^c) = (sine, 0,cos6') (5.4) 



in the right-hand-side of the SBE (5.2), we obtain com- 
plete cancellation. If we wish, we can follow HMH and 
separate the noise term into the contribution from feed- 
back (proportional to A) and that contribution present 
even without feedback (the rest). We interpret the latter 
stochasticity as being due to the quantum measurement 
we are making, with its underlying probabilistic nature. 
Obviously the fluctuation due to measurement is canceled 
by the feedback, as HMH point out. It is equally impor- 
tant that the deterministic dynamics are also canceled at 
this point. 

The story for the special case 6 — ■n/2 is quite different. 
For this case the SBEs are 
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(5.5) 



Here the three eigenvalues in Eq. ( 3.12| ) are clearly evi- 
dent. Both Zc and y,, will decay to zero (as required for 
= 7r/2), and their noise terms vanish at that point. By 
contrast, the equation for Xc is independent of the others, 
and is purely stochastic: 



dx,^ ^dW(t){l-xl). 



(5.6) 



Clearly the equatorial pure states with Xc = ±1 are 
stationary solutions to this problem. Also, the system 
will tend to one of these states. We can see this by cal- 
culating 



dF.[xl] = 7d^E [(1 - xlf] 



(5.7) 



which is always positive. That is, on average always 
increases. But it is also clear that Xc has no preference 
to go to either of these states. Hence they are not stable. 
The ensemble average x is unchanging under this evolu- 
tion. Thus a perturbation which moves the state from 
Xc = 1 to Xc = 1 — e say, will result in a proportion e/2 
of the states ending up at Xc = —1, and a proportion 
1 — e/2 ending up at Xc = 1 . 

We have illustrated these features by showing two typ- 
ical trajectories in Fig. ^(b). Once again, the initial evo- 
lution is highly erratic, but the system reaches a fixed 
point on a time scale of a few 7"^. However, with the 
same initial condition (the ground state), one trajectory 
ends up at Xc = 1 and the other at Xc = —1. 



C. Non-unit-efHciency 

In the unit-efficiency case the stationary solution of the 
master equation is (except for \9\ = 7r/2), a pure state. 
This is very special in that in means that every unravel- 
ing of the master equation as a SSE or SME must end up 
in this same pure state also. For non-unit-efficiency we 
have found the most pure stable state for each 9. In this 
case we must use a SME to unravel the master equation, 
since the conditioned state will not be pure in general, 
because of the inefficient detection. Since the determin- 
istic steady state is not pure (except for \9\ = tt), the 
quantum trajectories need not end up in this state. In- 
stead, the quantum state in an individual trajectory may 
continue to evolve stochastically even when the system is 
in steady state, and the equivalence to the deterministic 
evolution may hold only on average. On the other hand, 
it is also possible that the quantum trajectories do all 
end up in the deterministic steady state, since we expect 
the conditioned state to be mixed anyway. 
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It turns out that with the optimal values of a and A 
defined in Sec. IV, the actual behaviour is the first op- 
tion described above. That is, the system state continues 
to vary stochastically in the long time limit, but is con- 
strained so that the timc-averaged state is equal to the 
solution of the deterministic master equation. We show 
this in Bloch space Fig. ^ for 77 = 0.8 and = tt/6. We 
see that the amount of randomness in the system state 
in the long-time limit is quite large even for fairly high 
V- 

This result suggests another question: would a differ- 
ent choice for A be able to reduce, or even eliminate, the 
randomness in the steady-state quantum trajectory, even 
though it would necessarily be at the expense of the pu- 
rity of the deterministic stationary solution. [R ecal l that 
for a given A, a is still necessarily fixed by Eq. (4.9).] To 
test this idea we tried choosing A based not on maximiz- 



ing as in Eq. ( 4.10| ), but on minimizing 



NeiX) = 



777 



(5.8) 



That is, we minimize the noise terms in the SBE 
Eq. (5.2). Note that we have replaced the conditioned 
Bloch variables Xc etc. with the deterministic station- 
ary solutions Xss etc., and that the dependence of these 
stationary solutions on a and A add a further, implicit, 
dependence on A to iVe(A). This is a sensible procedure 
if the aim is realizable, and the noise in the solutions 
is reduced or eliminated so that the conditioned states 
are approximately or exactly equal to the deterministic 
stationary solution. 

It turns out that this procedure cannot significantly re- 
duce the amount of steady-state randomness in the quan- 
tum trajectories below that resulting from minimizing 
the deterministic stationary purity. In fact, for all values 
of rj we considered, the variation of A (as a function of 9) 
based on minimizing the noise was indistinguishable by 
eye from that based on maximizing the purity. This is 
not too surprising, but could not have been predicted a 
priori. 



VI. CONCLUSION 

We have given a rigorous analysis of the anti- 
decoherence feedback scheme proposed by Hofmann, 
Mahler and Hess They proposed modulating the 

driving of a two-level atom using the instantaneous ho- 
modyne photocurrent, in order to stabilize the atom in 
an arbitrary known pure state. We have shown that, for 
detection efficiency = 1, the pure states thus produced 
are stable. This is contrary to the conclusion of HMH, 
that only pure states in the lower half of the Bloch sphere 
would be stable. The one exception we found is for pure 



states on the equator. Although they are fixed points of 
the dynamics, they are not stable. A small perturbation 
away from one fixed pure state leads to a proportionally 
small fraction of the ensemble ending up in the diamet- 
rically opposite pure state. 

It is nevertheless possible to obtain an asymmetry be- 
tween the upper and lower halves of the Bloch sphere, 
reminiscent of the conclusion of HMH, if one considers 
detection efficiencies less than one. In this case, it is no 
longer possible to stabilize the system in a given pure 
state, so we choose the feedback and driving so that the 
solution of the master equation (including feedback) is 
as close as possible to a given pure state. We find that 
the purity (which measures this closeness) of states thus 
produced decays to zero as rj decreases to zero, for states 
in the upper half of the Bloch sphere. By contrast, those 
in the lower half do not decay to zero. This is readily 
understandable since in the absence of feedback (which 
is the situation which must prevail when the detection 
efficiency goes to zero) , the master equation with driving 
alone has stationary solutions in the lower- half plane with 
non-zero purity. The purity decays most rapidly with rj 
for states near the equator, which is unsurprising given 
the instability of states on the equator even for rj = 1. 

In the non-unit-efficiency case, the state of the system 
conditioned on the homodyne measurement results con- 
tinues to evolve stochastically even in the long-time limit, 
where the ensemble average evolution has reached the de- 
sired most-pure state. Moreover, it seems that any other 
choice of driving and feedback will result in more, not 
less, randomness in the steady-state quantum trajectory. 

Our results are significant in a number of ways. First, 
they show the power of the quantum trajectory and mas- 
ter equation techniques developed in Rcfs. [§j8[|^. Those 
techniques were particularly useful for illuminating sub- 
tle questions regarding the stability of pure states, and 
for treating inefficient detection. Second, the physical 
system (the two-level atom) may one day find applica- 
tion as a quantum bit in quantum information technol- 
ogy In that eventuality, the ability to stabilize the 
atom against dissipation in an arbitrary (known) pure 
state may be useful. Third, the system is a simple but 
non-trivial example of quantum feedback in a nonlinear 
system (the two- level atom). Thus the effectiveness of 
feedback, and in particular the influence of non-unit- ef- 
ficiency detection on this effectiveness, is of interest for 
what it may tell us about other more complicated non- 
linear systems. 

In this last context, it would be of interest to also con- 
sider the effect of non-Markovian feedback; that is, feed- 
back with a time delay or non-flat loop response function. 
This is much more difhcult to treat than Markovian feed- 
back because the Lindblad master equations derived in 
Refs. do not apply. Analytical solutions for non- 

Markovian feedback are possible for linear systems [^,^ • 
For a nonlinear system like the two-level atom, numeri- 
cal simulations, or novel analytical approaches, are nec- 
essary. This is an issue we plan to explore in future work. 
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